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As a stringent test of the nonrelativistic QCD factorization approach for inclusive 
heavy-quarkonium production in high-energy hadron collision, the polarized heavy- 
quarkonium production at the Tevatron is discussed. The polarization of prompt 
J/tp at the Tevatron is predicted within the nonrelativistic QCD factorization 
framework. The contribution from radiative decays of P-wave charmonium states 
decreases, but does not eliminate, the transverse polarization at large transverse 
momentum. The prediction agrees with measurements from the CDF collaboration 
at intermediate values of px, but disagrees at the large values of px measured. Re- 
cent preliminary data of heavy quarkonia from CDF and theoretical / experimental 
prospects related to the Run II heavy quarkonium physics are also discussed. 



As an effective field theory of QCD, the nonrelativistic QCD (NRQCD) 
factorization approachEl has had a great success in describing the inclusive pro- 
duction and decay of heavy quarkonium systems which could not be explained 
by a conventional phenomenological model, the so-called color-singlet model 
(CSM). The color-octet mechanism of NRQCD, which is absent in the CSM, 
presented a resonable solution EL to the anomalous production rate of ip' and 
J/ip at the Fermilab Tevatron El. In high-energy pp collisions, the dominant 
contribution ta the charmonium production rate at large px comes from gluon 
fragmentationu. At leading order in a s , a QQ pair with small relative momen- 
tum created by the virtual gluon is in a color-octet 3 Si state with the same 
transverse polarization as the almost on-shell gluon. The enhancement of the 
production cross section due to the gluon propagator makes the fragmenta- 
tion process dominate. Since the short-distance contribution has two powers 
of a s less than the CSM gluon fragmentation channel, the high-pr S-wave 
charmonium production rate should be dominated by the process from the 
transition of color-octet 3 Si state even though there is a suppression factor v A 
coming irom a double chromoelectric dipole transition in the non-perturbative 
processo. _ 

The NRQCD factorization approach to inclusive quarkonium productionlil 
also makes the remarkable prediction that in hadron collisions, J PC = 1~~ 
quarkonium states should be transversely polarized at large pr"- Due to the 
heavy-quark spin symmetry of the NRQCD interaction Lagrangian govern- 
ing the chromoelectric dipole transition and the fact that the transition of 
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the QQ pair into a physcal S— wave charmonium is dominated by a double 
chromoelectric dipole transition, the polarizarion of J/ip mesons produced at 
sufficiently large pt should be transversely polarized. Recent measurements at 
the Tevatmn by the CDF collaboration seem to be in contradiction with this 
predictionO. 

A convenient measure of the polarization is a = (T— 2L)/(T+2L), where T 
and L are the transverse and longitudinal polarization fractions in the hadron 
CM frame. The variable a describes the angular distribution of leptons from 
the decay of the J/tp with respect to the J/ip momentum. Beneke and Roth- 
stein studied the dominant-fragmentation mechanisms for producing longitudi- 
nally polarized 1 statesQ. For charmonium production at the Tevatron, one 
should also take into account the fusion contributions from parton processes 
ij —> cc + k since fragmentation does not yet dominate for most of the pt 
region. The polarization variable a for direct J /iL and direct ip' at the Teva- 
tron have been predicted by Beneke and Kramer 19 and by LeibovichB. They 
predicted that a should be small for pt below about 5 GeV, but then should 
begin to rise dramatically. On the contrary, the CDF data shows no sign of 
transverse polarization of direct ip' at large pr^. _ 

The CDF collaboration has also measured the polarization of prompt J/ipti 
(i.e. J/ip's that do not come from the decay of B hadrons). The advantage 
of this measurement is that the number of J/ip events is larger than for ip' 
by a factor of about 100. On the other hand, theoretical predictions of the 
polarization of prompt J/ip are complicated since the prompt signal includes 
J/ip^ that come from decays of the higher charmonium states Xcij Xc2, and 
ip 1 E3ca The polarization of J/tp from ip' not via \c is straightforward to 
calculate, since the spin is unchanged by the transition. The polarization of 
J/ip from Xc and of J/tp from ip' via \c is more complicated, because the Xcj's 
are produced in various spin states and they decay into J/ip through radiative 
transitions E3. Therefore, it is interesting to investigate the cascade effect in 
this problem. 

The NRQCD factorization formula for the differential cross section for 
the inclusive production of a charmonium state H of momentum P and spin 
quantum number A has the schematic form 

da H ^ = da c ^ (O^), (1) 

where the implied sum on n extends over all the color and angular momentum 
states of the cc pair. The cc cross sections do CCn , which are independent of 
H, can be calculated using perturbative QCD. All dependence on the state 
H is contained within the nonperturbative NRQCD matrix elements (ME's) 
(On X ^)- The Lorentz indices, which are suppressed in (|l|), are contracted 
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with those of da CCn to give a scalar cross section. The symmetries of NRQCD 
can be used to reduce the tensor ME's (On X ^) to scalar ME's (0%) that are 
independent of P and A. One may calculate the cross section for polarized 
quarkonium once the relevant scalar ME's are known. 

In pp collisions, the parton processes that dominate the cc cross section de- 
pend onpx- If Pt is of order m c , the fusion processes dominate. These include 
the parton processes ij — > cc + k, with i, j = g,q,q and q = u,d, s. At px much 
larger than m c , the parton cross sections are dominated by fragmentation pro- 
cesses. To get a consistent LO evaluation of polarization parameter a in this 
region, one should take care. Since we are interested in the LO prediction of 
a, we should first get a LO prediction of the longitudinal fraction L/(T + L). 
The LO contributions to L are order-o^ processes g — > cc-jf- g, while the only 
contribution to T is the order-a s process g ccg ( 3 Si)0. Inclusion of the 
longitudinal part in the fragmentation process is crucial in the calculation of 
a even though it has only a small contribution to the total cross section. The 
fragmentation function is evolved using the standard homogeneous timelike 
evolution equation. Since cc$( 3 Si) — > ip\(nS) is included in both the fragmen- 
tation and the fusion process, we interpolate between the fusion cross section 
at low pt and the fragmentation cross section at high p^"'"- In ah the other 
cc — > H channels, we use the fusion cross section. 

The color-singlet ME's can be determined _jihenomenologically from the 
decay rates for ip(nS) — > t + t~~ and \c2 — * 77 Ej. The color-octet ME's are 
phenomenologicaLparameters and they are determined from production data 
at the Tevatrontl'Ej. In-pase of Xa we have a good agreement with the CDF data 
on the rat.ia XcilXci More details on the analysis methods are explained 
in Refs.liatl There are numerous theoretical uncertainties in the polarization 
calculation. Ill our analysis, we allow the variations in PDF (CTEQ5L and 
MRST98LO) factorization and fragmentation scales (/*t/2 — charm 
quark mass (1.45 — 1.55 GeV), respectively, where \it = {^m^+p^) 1 / 2 . We also 
take into account the uncertainties in the ME's, (OsO'S'o)} and (Os( 3 Po)) 

Let us compare our results for a with the CDF data. We present our 
result in the form of an error band obtained by combining in quadrature all the 
theoretical errors described previously. In Fig. 1(a), we compare our result for 
direct ijj' as a function of px with the CDF dataB. But the error bars in the CDF 
data are too large to draw any definitive conclusions. We next consider the 
polarization variable a for prompt J /if}. The method fanxalculating cascade 
effect in prompt J/ip polarization is explained in Ref. 113. In Fig. 1(b), we 
compare our result for a as a function of pt with the CDF datacl. Our result 
for a is small around px = 5 GeV, but it increases with px- Our result is in 
good agreement with the CDF measurement at intermediate values of px, but 
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Figure 1: Polarization variable a vs. for (a) direct j/i' and (b) prompt J/?/) compared to 
CDF data. 



it disagrees in the highest pr bin, where the CDF measurement is consistent 
with 0. The solid lines in Fig. 1(b) are the central curves of a for direct J/ip 
and for J/ip from Xc- The a for direct J/tp is smaller than that for direct ip' , 
because the ME's are differentEj. In the moderate pt region, the contributions 
from ip' and from Xc add to give an increase in the transverse polarization of 
prompt J/ip compared to direct J/ip. In the high pt region, the contributions 
from ip' and Xc tend to cancel. 

The CDF measurement of the polarization of prompt J/tp presents a seri- 
ous challenge to the NRQCD factorization formalism for inclusive quarkonium 
production. The qualitative prediction that a should increase at large pt seems 
inescapable. However, it is still worthwhile to investigate the NLO effect in the 
fragmentation processes, since the dramatic discrepancy appears at large px- 
Recently, Braaten and Lee performed the full NLO calculation of the calor- 
octet 3 Si gluon fragmentation function for polarized heavy quarkonium 113. It 
is not sufficient to get NLO prediction of a for prompt J/ip, but one can first 
study the production rate at large pr in NLO accuracy. Full NLO calculation 
of fusion process in NRQCD framework is also necessary in order to get more 
precise values for the ME's. 

In Run II of the Tevatron, the data sample of heavy quarkonium should be 
more than an order of magnitude larger than in Run I. This will allow for both 
the production rate and the polarization to be measured with higher precision 
and out to larger values of pr- The progress of the experimental analysis in 
Run I will advance further. The NRQCD ME's for both the charmonium and 
bottomonium systems will be determined with higher accuracy. Especially, the 
accuracy of the color-octet 3 Si ME for the ip' meson will be greatly increased 
since one can probe the high-p-r region where fragmentation dominates and 
which currently has little statistics. The measurements of the polarization 
parameter a for prompt J/ip and ip' will be more precise. It will also be 
possible to measure a's for direct J/ip and J/ip from Xcj since the experimental 
analysis technique to separate the prompt signal is well developed already. The 
polarization parameter for J/ip from Xcj decay, which is predicted in Fig. 1(b) 
as well as that for direct J/ip can also be measured. Since CDF successfully 
separated the Xci and Xc2 signal in their total cross section measurement, the 
Pt distributions of both Xci and Xci can be measured separately. The heavy- 
quark spin symmetry which has been used for the simplification of Xcj ME's 
can be tested with higher accuracy. 

If the polarization data continues to disagree with the predictions of the 
NRQCD factorization approach, it would indicate a serious flaw in our un- 
derstanding of inclusive charmonium production. If it is due to the fact that 
m c is too small, it might be interesting to test the polarization of T. The 
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recent measurement of the T production rate by CDF will allow for a better 
understanding of the NRQCD matrix elements for the T system without many 
assumptions which were necessary due to the lack of statistics. The theoretical 
prediction of the T polarization is needed both for Run I and Run II. We might 
also have a chance to probe the gluon fragmentation dominance in high-p^ Y 
production in Run II even though the mass of T is three times larger than that 
of J/ip. 

We expect that Run II of the Tevatron will definitely produce significant 
information which will help us to understand charmonium and bottomonium. 
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